Few-Body Systems 0, 1-7 (2008) p 

bystems 



© by Springer- Vcrlag 2008 
Printed in Austria 



in 
o 
o 



On 



in 
o 



How To Classify 3-Body Forces — And Why 



Harald W. GrieBhammer* ** 



> 

q ' Institut fiir Theoretische Physik III, Universitat Erlangen, D-91058 Erlangen, Germany 

and Institut fiir Theoretische Physik (T39), TU Miinchen, D-85747 Garching, Germany 



Abstract. For systems with only short-range forces and shallow 2-body 
bound states, the typical strength of any 3-body force in all partial-waves, 
including external currents, is systematically estimated by renormalisation- 
. group arguments in the Effective Field Theory of Point-Like Interactions. The 

underlying principle and some consequences in particular in Nuclear Physics 
are discussed. Details and a better bibliography in Ref. [1]. 



1 Introduction 



O 

^ | The Effective Field Theory (EFT) of Point-Like Interactions is a model- 

independent approach to systems without infinite-range forces in Atomic, Molec- 
ular and Nuclear Physics at very-low energies with shallow real or virtual 2-body 
bound-states ("dimers"), see e.g. [2,3] for reviews. When the size or scattering 
' length a of a 2-body system is much larger than the size (or interaction range) 

R of the constituents, a small, dimension-less parameter Q = -§ allows to clas- 
sify the typical size of neglected corrections at nth order beyond leading order 
(N n LO) as about Q n . For example, a ~ 104 A and R ~ 1 A in the 4 He2 molecule, 
i.e. Q « yqq, while a 4.5 fm and R 1.5 in the deuteron, i.e. Q « | in the 
"pion-less" EFT, EFT(^f), where pion-exchange between nucleons is not resolved 
as non-local. Thus, the detailed dynamics on the "high-energy" scale R can vary 
largely: For example, attractive van-der-Waals forces oc -\ balance in 4 He2 a 
repulsive core generated by QED; but in Nuclear Physics, one-pion exchange 
oc 3 ] is balanced by a short-range repulsion whose origin in QCD is not yet 
understood. It is a pivotal advantage of an EFT that it allows predictions of 
pre-determined accuracy without such detailed understanding - as long as one 
is interested in low-energy processes, i.e. Physics at the scale a, and not R. Even 
when possible (as in QED - at least at scales > 1 fm), EFTs reduce numerically 
often highly involved computations of short-distance contributions to low-energy 
observables by encoding them into a few simple, model-independent constants of 
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contact interactions between the constituents. These in turn can be determined 
by simpler simulations of the underlying theory, or - when they are as in QCD 
not (yet) tenable - by fit to data. Universal aspects of few-body systems with 
shallow bound-states are manifest in EFTs, with deviations systematically cal- 
culable. In 2-body scattering for example, the EFT of Point-Like Interactions 
reproduces the Effective- Range Expansion, but goes beyond it in the systematic, 
gauge-invariant inclusion of external currents, relativistic effects, etc. 

Take 3-body forces (3BFs): They parameterise interactions on scales much 
smaller than what can be resolved by 2-body interactions, i.e. in which 3 particles 
sit on top of each other in a volume smaller than R 3 . Traditionally, they were 
often introduced a posteriori to cure discrepancies between experiment and the- 
ory, but such an approach is of course untenable when data are scarce or 1- and 
2-body properties should be extracted from 3-body data. But how important are 
3BFs in observables? The classification in EFTs rests on the tenet that a 3BF 
is included if and only if necessary to cancel cut-off dependences in low-energy 
observables. I outline this philosophy and its results in the following, finding that 
- independent of the underlying mechanism - 3BFs behave very much alike in 
such disparate systems as molecular trimers and 3-nucleon systems, but do not 
follow simplistic expectations. 

2 Construction 

In the Faddeev equation of particle-dimer scattering without 3BFs, Fig.^ the S- 
wave 2-body scattering amplitude is given by the LO-term of the Effective- Range 
Expansion. This dimer and the remaining particle "interact" via /Q, the one- 
particle propagator projected onto relative angular momentum I. Even for small 




N n LO 



Figure 1. Left: integral equation of particle-dimer scattering. Right: generic loop correction 
(rectangle) at N n LO. Thick line (D): 2-body propagator; thin line (/G): propagator of the 
exchanged particle; ellipse: LO half off-shell amplitude t\ (p). 

relative on-shell momenta k between dimer and particle, we need the scattering 
amplitude t^ip) for all off-shell momenta p to determine its value at the on- 
shell point p = k, and hence in particular for p beyond the scale on which a 
description in terms of point-like constituents is tenable. It is therefore natural to 
demand that all low-energy observables on a scale k ~ - are insensitive to derails 
of the amplitude at p > j^, namely to form and value of the regulator, form- 
factor or cut-off chosen. If not, a 3BF must soak up the dependence. In an EFT, 
this is the fundamental tenet: Include a 3BF if and only if it is needed as counter- 
term to cancel divergences which can not be absorbed by renormalising 2-body 
interactions. Thus, only combinations of 2- and 3BFs are physically meaningful. 
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With the cut-off variation of the 3BF thus fixed, the initial condition leads to 
one free parameter fixed from a 3-body datum or knowledge of the underlying 
physics. 3BFs are thus not added out of phenomenological needs but to guarantee 
that observables are insensitive to off-shell effects. 

A Mellin transformation t^(p) oc p^M^ 1 solves the equation for p^> k,^. 
The spin-content is then encoded only in the homogeneous term: A = — | for 
3 nucleons with total spin |, or for the totally spin and iso-spin (Wigner-) anti- 
symmetric part of the spin-^-channel; A = 1 for 3 identical spin-less bosons and 
the totally spin and iso-spin (Wigner-) symmetric part of the spin-i-channel. The 
asymptotic exponent s;(A) has to fulfil Re[s] > —1, Re[s] ^ Re[l + 2], and 



1 = (-1) 



2l _, A r p+^+l] r r. 



2Z+31 



2-^1 



l + s + 1 l-s + l 21 + 3 1 
2 ' 2 ' 2 ' 4 



This result was first derived in the hyper-spherical approach by Gasaneo and 
Macek [4] 1 . The asymptotics depends thus only but crucially on A and I. Relevant 
in the UV-limit are the solutions for which Re[s + 1] is minimal. 

At first glance, we would expect the asymptotics to be given by the asymp- 
totics of the inhomogeneous (driving) term: t^\p) 6c -Jr^, i.e. si(X) =1 + 1. 
However, we must sum an infinite number of graphs already at leading order. As 
Fig- ID shows, this modifies the asymptotics considerably. 




-1 01234 -1 01234 

1 1 



Figure 2. The first two solutions s;(A) at A = 1 (left) and A = — i. Solid (dotted): real 
(imaginary) part; dashed: simplistic estimate. Dark/light: first/second solution. Limit cycle 
and Efimov effect occur only when the solid line lies below the dashed one, and Im[s] ^ 0. 



How sensitive are higher-order corrections to the UV-behaviour of vJ (p)? The 
2-body scattering-amplitude is systematically improved by including the effective 
range, higher partial waves etc. Corrections to 3-body observables (including 
partial- wave mixing) are found by perturbing around the LO solution as in Fig.Q 
Most sensitive to unphysically high momenta is each correction at N ra LO which 
is proportional to the nth power of loop momenta. The question when it becomes 
cut-off sensitive is now rephrased as: When does the correction diverge as the 
cut-off is removed, i.e. when is its superficial degree of divergence non-negative? 

1 My apologies to the authors that I found this reference only after [1] was published. 
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The answer by simply counting loop momenta in the diagram: 



Re[n- S/ (A)-^(A')] >0 . 




We therefore find at which order the first 3BF is needed just by determining when 
a correction to the 3-body amplitude with only 2-body interactions becomes 
dependent on unphysical short-distance behaviour. 

It is instructive to re-visit these findings in position space. The Schrodinger 
equation for the wave-function in the hyper-radial dimer-particle distance r, 



looks like the one for a free particle with centrifugal barrier. One would thus 
expect si = I + 1 (hyper-spherical co-ordinates!). It had however already been 
recognised by Minlos and Faddeev that the centrifugal term is for three bosons 
(A = 1) despite expectations attractive, so that the wave-function collapses to 
the origin and seems infinitely sensitive to very-short-distance physics. In order 
to stabilise the system against collapse - or, equivalently, remove dependence on 
details of the cut-off -, a 3BF must be added - or, equivalently, a self-adjoint ex- 
tension be specified at the origin i.e. a boundary condition for the wave-function 
must be fixed by a 3-body datum. On the other hand, 3BFs are demoted if 
si > I + 1: The centrifugal barrier provides more repulsion than expected, and 
hence the wave-function is pushed further out, i.e. less sensitive to details at dis- 
tances r < R where the constituents are resolved as extended. Birse confirmed 
this recently by a renormalisation-group analysis in position-space [5]. 

3 Consequences 

About half of the 3BFs for I < 2 are weaker, half stronger than one would expect 
simplistically, see Table ^ The higher partial-waves follow expectation, as the 
Faddeev equation is then saturated by the Born approximation. The S-wave 3BF 
of spin-less bosons is stronger, while the P-wave 3BF is weaker. 

That the first S-wave 3BF appears already at LO leads to a new 
renormalisation-group phenomenon, the "limit-cycle". It explains the Efimov 
and Thomas effects, and universal correlations e.g. between particle-dimer scat- 
tering length and trimer binding energy (the Phillips line). In general, it appears 
whenever the kernel of the integral equation not compact, i.e. Im[s] ^ and 
|Re[s]| < Re[^ + 1]. We finally note that the power-counting requires a new, inde- 
pendent 3BF with 21 derivatives to enter at N^LO and provides high-accuracy 
phase-shifts in atom-dimer and nucleon-deuteron scattering, and loss rates close 
to Feshbach resonances in Bose-Einstein condensates, see e.g. [2, 3, 6] for details. 

Demotion might seem an academic dis-advantage - to include some higher- 
order corrections which are not accompanied by new divergences does not im- 
prove the accuracy of the calculation; one only appears to have worked harder 
than necessary. However, demotion is pivotal when one wants to predict the ex- 
perimental precision necessary to dis-entangle 3BFs in observables, and here the 
error-estimate of EFTs is crucial. In 4 He-atom-dimer scattering, where Q ~ jirr, 




(3) 
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bosons fcrmions 


naive dim. analysis 
Refs/fA) + si'(X')] 


simplistic 
l + l' + 2 




typ. size 
if n n ~ — 


S~ S S~ S 
2 S- 4 D 


LO 

N 31 LO 


N 2 LO 
N 4 LO 


prom, 
prom. 


100% (10%) 
3% (1%) 


P_P ^p_^p 

2p_2p 4p_4p 
2p_4p 


N bY LO 
N 35 LO 
N 4 - 6 LO 


N 4 LO 


dem. 
prom, 
dem. 


0.2% 

2% (1%) 
0.6% 


4 S- 4 S 
4 S- 2 D 
4 S- 4 D 


N 4 - 3+2 L0 
N 50 LO 
N 53 LO 


N 2+2 LO 
N 4 LO 


dem. 
dem. 
dem. 


0.1% 

0.4% (1%) 
0.3% 


higher 


~ as simplistic 


N i+r+2 LO 







Table 1. Order of the leading 3BF, indicating if actual values are stronger ("prom.") or 

weaker ("dem.") than the simplistic estimate. Last column: typical size of 3BF in EFT(^f); in 
parentheses size from the simplistic estimate. 

only high-precision experiments can however reveal contributions from 3BFs be- 
yond the one found in the S-wave. 

On the other hand, Q ~ | in EFT(ff) of Nuclear Physics. Now, the demotion 
or promotion of 3BFs makes all the difference whether an experiment to deter- 
mine 3BF effects is feasible at all. For example, the quartet-S wave scattering- 
length in the neutron-deuteron system sets at present the experimental uncer- 
tainty in an indirect determination of the doublet scattering length, which in 
turn is well-known to be sensitive to 3BFs. Its value in EFT(^) at N 2 LO, 

a( 4 S|) = [5.09(LO) + 1.32(NLO) - 0.06(N 2 LO)] fm = [6.35 ± 0.02] fm , (4) 

converges nicely and agrees very well with experiment, [6.35 ± 0.02] fm. The 
theoretical accuracy by neglecting higher-order terms is here estimated conser- 
vatively by Q ~ | of the difference between the NLO- and N 2 LO-result. Tabled 
predicts that the first 3BF enters not earlier than N 6 LO. Indeed, if the theo- 
retical uncertainty continues to decreases steadily as from NLO to N 2 LO, an 
accuracy of ±(^) 3 X 0.02 fm < ±0.001 fm with only 2-nucleon scattering data as 
input can be reached in calculations. This is comparable to the range over which 
modern high-precision potential-model calculations differ: [6.344 . . . 6.347] fm. If 
the 3BF would occur at N 4 LO as simplistically expected, the error by 3BFs 
would be (|) 4 x 0.02 fm ~ 0.007 fm, considerably larger than the spread in 
the potential-model predictions. Differential cross-sections and partial-waves are 
also in excellent agreement with much more elaborate state-of-the-art potential 
model calculations at energies up to 15 MeV, see e.g. [6]. 

The cross-section of triton radiative capture nd — > t*y at thermal energies 
provides another example [7]. Nuclear Models give a spread of [0.49 . . . 0.66] mb, 
depending on the 2-nucleon potential, and how the Z\(1232) as first nucleonic 
excitation is included. On the other hand, a process at 0.0253 eV [sic] incident 
neutron energy and less than 7 MeV photon energy should be insensitive to de- 
tails of the deuteron wave-function and of a resonance with an excitation energy 
of 300 MeV. Indeed, the power-counting of 3BFs applies equally with external 
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currents, only that the higher-order interaction in Fig. ^ includes now also the 
momentum- or energy-transfer from the external source as additional low-energy 
scales. As no new 3BFs are needed up to N 2 LO to render cut-off independent 
results, the result is completely determined by simple 2-body observables as 

a tot = [0.485(LO) + 0.011(NLO) + 0.007(N 2 LO)] mb = [0.503 ± 0.003] mb , (5) 

which converges and compares well with the measured value, [0.509 ± 0.015] mb. 
The cross-section relevant for big-bang nucleo-synthesis (E n ~ 0.020 . . . 0.4 MeV) 
is also in excellent agreement with data [8]. 

4 Conclusions 

With these findings, the EFT of 3-body systems with only contact interactions 
is a self-consistent, systematic field theory which contains the minimal number 
of interactions at each order to render the theory renormalisable. Each 3-body 
counter-term gives rise to one subtraction-constant which is fixed by a 3-body 
datum. Table ^ sorts the 3BFs by their strengthes, their symmetries and the 
channels in which they contribute at the necessary level of accuracy. Amongst the 
host of applications in Nuclear Physics are triton and 3 He properties, reactions in 
big-bang nucleo-synthesis, neutrino astro-physics, the famed nuclear A^-problem, 
and the experimental determination of fundamental neutron properties. 

The method presented here is applicable to any EFT in which an infinite 
number of diagram needs to be summed at LO, e.g. because of shallow bound- 
states. One example is Chiral EFT, the EFT of pion-nucleon interactions. Only 
those local A-body forces are added at each order which are necessary as counter- 
terms to cancel divergences at short distances. This mandates a careful look at 
the ultraviolet-behaviour of the leading-order, non-perturbative scattering am- 
plitude. It leads at each order and to the prescribed level of accuracy to a cut-off 
independent theory with the smallest number of experimental input-parameters. 
The power-counting is thus not constructed by educated guesswork but by rig- 
orous investigations of the renormalisation-group properties of couplings and 
observables using the methodology of EFT. 
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